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ABSTRACT 

Context. Ground based Cherenkov telescope systems measure astrophysical y-ray emission against a background of cosmic-ray induced air 

showers. The subtraction of this background is a major challenge for the extraction of spectra and morphology of y-ray sources. 

Aims. The unprecedented sensitivity of the new generation of ground based very-high-energy y-ray experiments such as H.E.S.S. has lead 

to the discovery of many previously unknown extended sources. The analysis of such sources requires a range of different background 

modelling techniques. Here we describe some of the techniques that have been applied to data from the H.E.S.S. instrument and compare their 

performance. 

Methods. Each background model is introduced and discussed in terms of suitability for image generation or spectral analysis and possible 
caveats are mentioned. 

Results. We show that there is not a single multi-purpose model, different models are appropriate for different tasks. To keep systematic 
uncertainties under control it is important to apply several models to the same data set and compare the results. 
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1. Introduction 

Ground based very-high -energy y-ra y telescop e sys- 
tems such as H. E.S.S. faintonl l2004 . MAGIC JLorenzl 
l2004. VERITAS IWeekes et all bOodT and CANGAROO- 
3 IlKubo et alJ I2004I) have greatly increased the sensitivity 
of the Atmospheric Cherenkov technique. However, these 
instruments can only reach their full potential if systematic 
effects are brought fully under control. A major challenge for 
experiments of this type is the subtraction of the background 
of non-y-ray induced air showers. This background can be 
dramatically reduced using image-shape selection criteria, but 
cannot be removed completely. The background above a few 
hundred GeV is dominated by hadronic cosmic-ray showers, 
with cosmic-ray electrons increasingly important at low 
energies and after tight image selection cuts. The subtraction 
of this background is the main source of systematic errors and, 
if not done correctly, can even produce an artificial source. 

For single telescope instruments (for exa mple the pioneer- 
ing Whipple telescope IWeekes et alJll989l) ') the classical ap- 
proach to background subtraction was the ON/OFF observing 
mode. In this mode observations (runs) centred on the target 
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source are interspersed with equal-length observations of an 
empty field at equal zenith angle (typically a region offset in 
Right Ascension by 30 minutes). The background is assumed 
to be equal in the two runs, the difference between them pro- 
vides a measurement of the y-ray signal. A major drawback 
of this approach is that only half of the available dark time is 
spent ON-sou rce. The wobble-m ode pioneered by the HEGRA 
collaboration JPaum et alJll997t) avoids this problem by keep- 
ing the targeted source region in the field of view (FoV) at all 
times, with an alternating offset relative to the system's point- 
ing direction (for point sources typically ±0.5° in Declination). 
A background estimate (OFF data) for the source region (ON 
data) can then be derived from a region on the opposite side of 
the FoV from the same run as the ON data. 

For wide FoV instruments the probability of serendip- 
itous detection of non-targeted sources rises dramatically 
(particularly for observations close to the Galactic plane). 
The discovery of such sou rces (for example HESS J 1303- 
631 JAharonian et al]l2()0 5a)') demonstrates the need for back- 
ground models that provide background estimates for the 
whole FoV of the instrument. Any systematic survey of a 
whole sky region also requires such models. The importance 
of surveys has been demonstrated by the recent discovery 
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Fig. 1. The variation of the radial system acceptance function with zenith angle for two different sets of cuts. The distributions 
show the squared angular distance ifr 2 between reconstructed event directions and pointing direction of the telescope system 
generated from OFF-somce data. The curves are generated by smoothing the one-dimensional acceptance histograms. They are 
arbitrarily normalised to 1 at the system pointing direction, corresponding to ifr 2 - 0. Left: The standard loose set of cuts (labelled 
std) is shown. It is mainly used for the determination of y-ray-ray spectra. Right: The hard cuts employing a larger cut on the 
minimum image amplitude are shown. This configuration is typically used for morphology studies. 



of more than ten new y-ra y sources in a survey o f the in - 
ner Galaxy with H.E.S.S. jAharonian et alJ l2005bl l2006al) . 
More over, extended y-ray s ources such as RX J 17 13.7- 
3946 JAharonian et al.ll2004 | 2006bl) and RX J0852.0-4622 
(Vela Junior) ( Aha ronian et al.ll200 5c) present additional diffi- 
culties for background subtraction. They require correct back- 
ground modelling over a region of the sky that is substantially 
larger than the source itself. 

Many different approaches to background modelling are 
possible and have been applied in the analysis of data from the 
H.E.S.S. instrument. These models have different strengths and 
weaknesses and it is usually desirable to apply several of them 
in the analysis of any y-ray source to get a handle on the sys- 
tematic uncertainties connected to the background determina- 
tion. In this paper we aim to describe some of these approaches. 
While we will be using H.E.S.S. data as the basis for our exam- 
ples it should be noted that the techniques described here apply 
generally to VHE y-ray instruments. Where necessary we will 
point out features that are specific to the analysis applied or to 
H.E.S.S. and its large FoV. 

An outline of this paper is as follows: We begin in section|2] 
with a general introduction of background modelling and in- 
vestigate in detail the properties of the system acceptance, as 
it is a key starting point for almost all background models. We 
then move on to a description of individual background mod- 
els. In section|5]we compare the results of different approaches 
applied to particular data examples. In section^] the impact of 
bright stars on observations using atmospheric Cherenkov tele- 
scopes is demonstrated. Section|5]summarises the strengths and 
weaknesses of the different background models and describes 
to which purpose they are best suited. 



2. Background Modelling 

With cuts on image-shape parameters the cosmic-ray back- 
ground can be reduced by a factor of ~100, resulting in a 
signal-to-background ratio for a strong point source lik e the 
Crab nebula on the order of 1:1 ( Aharonia n et al. ll200 6cl The 
remaining background of y-ray-like events must be estimated 
to derive the significance of any possible y-ray signal. Given 
a number of counts N on in a test region, and N s counts in a 
background control region, the y-ray excess is defined as 

N exceS s=N OD -aN oS . (1) 

The parameter a is a normalisation factor which accounts for 
solid-angle, exposure-time, zenith-angle, and acceptance dif- 
ferences between the test region and the background control 
region. In general, it is the ratio of the effective (acceptance- 
weighted) exposure integrated in time and angular space over 
the signal and background region, usually referred to as the ON 
and OFF regions. It can generally be defined as: 

L AlJfi*, 6 r , <p z , t) d0 x d0 y dcf> z dt 
Iff A off(^' y> fa' f ) de * d0 y d & dt ' 

A y on off is the system acceptance of y-ray like events and depends 
on the position (6 X , y ) in the FoV and the zenith angle <f> z of 
observations. Additionally, different exposure times t for ON 
and OFF region have to be taken into account. Note that if ON 
and OFF region are different in size and shape, the integration 
of the system acceptance in Eq. 13 must take this into account. 
Given a number of ON and OFF counts and a, the statistical 
significance (S) o f the excess is ty pically calculated following 
the prescription of lLi & Ma ( 1983, equation 17). 
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The task of a background model is to provide the quantities 
A^ ff and a. A choice of background regions such that a <K 1 
(achievable, e.g., by choosing much larger OFF than ON re- 
gions) results generally in higher statistical significance, be- 
cause background fluctuations are reduced, but may also result 
in increased systematic errors. The principle difficulty in deriv- 
ing a background estimate is the determination of the correct 
value of a. Since proper control over the system acceptance is 
crucial for this purpose, we investigate below the acceptance of 
H.E.S.S.. 

2.1. Cosmic-Ray System Acceptance after y-ray Cuts 

The system acceptance is defined as the probability of accept- 
ing, after triggering, analysis cuts and y-ray selection, a back- 
ground event reconstructed at a certain position in the system 
FoV and with a certain energy. For most background models 
some knowledge of the system acceptance is required to gen- 
erate an image of y-ray excess events or calculate significances 
of arbitrary positions in the FoV. In general, the acceptance de- 
pends on: 

□ the position in the FoV, particularly the distance to the op- 
tical axis 

□ the zenith (<f> z ) and azimuth §. dz angle of observations, due 
to the influence of the Earth's magnetic field on the shower 
development in the atmosphere and the rotational asymme- 
try of the telescope system. 

□ the reconstructed primary energy 

□ the time of observation, due to possible changes in the sys- 
tem configuration and aging of mirrors 

□ the sky coordinates viewed, due to the night-sky back- 
ground light level. 

In most cases it is a reasonable assumption that the accep- 
tance is radially symmetric (the validity of this assumption 
is discussed later). It is generated in a one-dimensional fash- 
ion as the number of background events as a function of the 
(squared) angular distance between reconstructed event direc- 
tion and system pointing direction. It can either be determined 
on a run-by-run basis from the data set under analysis or be 
extracted from observations without significant y-ray emission 
in the FoV (OFF runs). In the latter case it is assumed that the 
system acceptance is identical for the ON and OFF runs. In the 
former case one may face two problems, y-ray contamination 
by a source, and lack of statistics. For a typical data run lasting 
28 minutes, recorded at moderate zenith angles, the number of 
events after y-ray selection cuts (available for the determination 
of the acceptance) is (9(10 4 ) with cuts for spectral analysis, and 
as low as O(l0 3 ) with cuts for morphology studies (see below 
for a description of the analysis cuts). 

Here we use 220 hours of H.E.S.S. observations without 
significant y-ray sources in the FoV to obtain a model of the ra- 
dial system acceptance. These reference observations are sub- 
divided into zenith-angle bands. Events passing y-ray cuts (i.e. 
y-ray-like background events) are then binned according to the 
squared angular distance between the reconstructed event di- 
rection and the system's pointing direction. Figure^iHustrates 



the dependence of the acceptance on the zenith angle of obser- 
vations and analysis cuts. The shallow central peak and rapid 
decline towards larger distances, stems from the analysis based 
on image (Hillas) parameters, where a cut on the distance be- 
tween image centre-of-gravity and the camera centre is applied 
to avoid truncation effects at the camera edge. Due to the finite 
camera size, edge effects are inherent and will always appear 
independent of the exact analysis applied. 

As can be seen from Fig.[Q 2° away from the system cen- 
tre, the y-ray acceptance at moderate zenith angles decreases 
to 20% - 50% of the peak value, depending on analysis cuts. 
In addition a smooth variation with zenith angle is apparent. 
With increasing zenith angle, the system acceptance broadens, 
an increasing fraction of events with directions further away 
from the system pointing direction is detected. This is a direct 
consequence of the fact that with increasing zenith angle the 
shower maximum is increasingly further away from the tele- 
scope system causing a broadening of the Cherenkov light-pool 
on ground and hence an enlarged phase space for events with 
large inclination angles. When comparing the average curve for 
any given zenith-angle band to the radial acceptance in differ- 
ent fields of view, observed at the same altitude, the scatter is 
relatively small, less than 3% within 1° of the observation po- 
sition and less than 10% out to 3°. It is therefore justifiable to 
use OFF data taken in different fields of view to determine a 
model of the system acceptance. 

The influence of analysis cuts is also apparent in Fig. ^ 
The two sets of cuts used throughout this paper are labelled 
std and hard. The first set includes a cut on the minimum am- 
plitude of each camera image at 80 photo-electrons (p.e.) and 
is optimised for the determination of source spectra. The sec- 
ond set uses a cut on the image amplitude at 200 p.e., and pro- 
vides better background suppression and superior angular res- 
olution. It is therefore normally used for source searches and 
image generation (more detailed descripti ons of the H.E.S.S . 
anal ysis techniques may be found in lAharonian et all (|2£)05d) 
and lAharonian et all J2006bl) '). The larger cut on the minimum 
image size results in curves which exhibit a generally less pro- 
nounced peak and a less rapid decline towards large distances. 
There is an increased fraction of events with large inclination 
angles with respect to the system pointing direction. 

The azimuth dependence of the radial system acceptance 
is small and therefore neglected here: when sub-dividing data 
taken in a narrow zenith-angle band into azimuth bins (say 
North, East, South, and West), only marginal differences occur 
at the few-percent level. The energy dependence of the accep- 
tance is much stronger, greatly complicating the use of back- 
ground models that require an acceptance correction for spec- 
tral analysis. This is illustrated in Fig.|2](left) where the energy 
dependence for a zenith angle range from 0° to 20° is plotted. 
The curves shown correspond to three different energy bands, 
E < 0.6 TeV , 0.6 TeV < E < 1.4 TeV, and 1.4 TeV < E. For 
relatively small energies the acceptance declines rapidly with 
increasing offset. For large energies the shape is completely 
different. High-energy showers result in large Cherenkov light- 
pool radii on ground. Therefore, as already mentioned, more 
events with large angular offsets start to trigger the array. In ef- 
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Fig. 2. Left: The energy dependence of the system acceptance is demonstrated. For the three different energy bands shown, the 
shape of the acceptance broadens dramatically with increasing energy. Right: Plot to test for a zenith-angle dependent FoV 
gradient (note the change of scales compared to the left-hand side). For OFF data taken at zenith angles of 40° to 45°, the 
symmetry of the system acceptance along an axis running parallel to the zenith direction is investigated. For this purpose we fill 
all events into a histogram in bins of zenith-angle difference A(p z between the event direction and the system pointing direction. 
To test the symmetry with respect to the pointing direction, the ratio is calculated from the number of events in bins on either 
side of the central A0 Z = bin, equally far away from this bin. In the plot we show this ratio versus \A<p z \. The data are split 
according to the azimuth of the observations: to the north (black crosses) and to the south (green crosses). Note that all corrections 
accounting for mis-pointing of the telescopes and for zenith-angle dependent trigger rates (see main text) are included here. 



feet, for energies beyond 1 .4 TeV, the acceptance is almost flat 
out to a distance of 2° from the system pointing direction. 

As previously mentioned, in most cases the system accep- 
tance is assumed to be radially symmetric. The most intuitive 
cause of deviations from radial symmetry is a zenith-angle de- 
pendent linear gradient across the FoV. The larger the zenith 
angle of observation, the larger the effective energy threshold of 
the system, due to the increasing absorption of showers. Since 
the energy spectrum of the cosmic -ray background is rather 
steep, the trigger rate, and thus the event rate, of the system 
decre ases smoothly with increasing zenith angle JFunk et all 
l2004 . Hence, in the H.E.S.S. FoV of « 5° significant variations 
of the system acceptance along the zenith axis may occur, and 
indeed such variations are observed. Depending on zenith an- 
gle the peak system acceptance extracted from OFF runs does 
not coincide with the nominal centre of the system pointing di- 
rection. It is found to be slightly shifted towards smaller zenith 
angles, away from the pointing direction. To account for this 
effect, the zenith-angle dependence of the shift is determined: 
for small zenith angles (~ 10°) the shift is negligible (< 0.01°), 
at moderate 30° it is on the order of 0.03° and exceeds 0.13° 
for very large zenith angles beyond 60°. A parametrisation is 
used to correct the nominal centre of the system pointing di- 
rection in each OFF run. Figure |2] (right) explores remaining 
deviations from radial symmetry along the zenith axis. For this 
purpose, OFF runs are processed in zenith-angle bands, stor- 
ing event distributions as a function of the zenith-angle differ- 
ence between the pointing direction of the system and the event 
direction A<p z . From the resulting distributions (ranging from 
A<p z = ±3.5°) a ratio is created dividing the number of events in 
a certain bin on the positive A0 z -side by the number of events in 



the corresponding bin on the negative side to test for symmetry 
around the zenith angle of observation. The resulting distribu- 
tion is shown in Fig. (right) for two azimuthal bands (north 
and south) to search for effects related to the Earth's magnetic 
field. If there was no zenith-angle dependence, the ratio would 
be 1 for the whole FoV. There seems to be a residual distortion 
of the system acceptance along the zenith axis, in the direc- 
tion one would expect from the trigger-rate variation. Larger 
zenith angles have smaller event numbers, indicating a slight 
under-correction of this effect. However, remaining deviations 
are estimated to be less than 5%. Within statistics, there is no 
North-South effect apparent, the event-ratio distributions are in 
reasonable agreement with each other. 

The validity of the simplifying assumption that the sys- 
tem acceptance is radially symmetric can be verified with 
data. For that purpose such a symmetric model acceptance 
derived from OFF data can be compared to the acceptance 
of a single data set without y-ray source, for example the 
H.E.S.S. data from the 2003-20 04 observation campaign of 
SN 1006 jAharonian et alJ2005eh (which shows no evidence of 
y-ray emission). Reconstructed directions of y-ray-like events 
are plotted in a coordinate system centred on the system point- 
ing direction in the "altitude-azimuth" system (the system is 
also referred to as nominal system). Accumulating events from 
different data runs one obtains in this way a sample two- 
dimensional y-ray acceptance map. This can be compared to 
a model acceptance map derived by choosing, for each run, 
the one-dimensional radial acceptance for the corresponding 
zenith angle (cf. Fig.Q. The acceptance is then rotated in the 
nominal system and accumulated for all runs yielding an accep- 
tance model which can be compared to the system acceptance 
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Fig. 3. Plots of the system acceptance are shown for observations of the supernova remnant SN 1006, compared to the radially 
symmetric model acceptance determined from OFF runs. Upper panel: Plots labelled "Data" and "Model" are the (arbitrarily 
normalised) acceptances determined in the nominal system on a run-by-run basis. Positive x-direction corresponds to positive 
altitude, positive y-direction to positive azimuth. The parallel horizontal and vertical lines define bands (0.8° wide) used to 
produce slices for a one-dimensional comparison, shown in the lower panel: The two plots show projections along x and y 
through data (red crosses) and model map (black lines) within the thick bands indicated in the upper panel. 



deduced from the data set. The result, derived from 6.3 hours 
(after dead-time correction) of (4-telescope) H.E.S.S. observa- 
tions of SN 1006, is shown in Fig. [3] There is general agree- 
ment between data and model acceptance. Remaining differ- 
ences are at the few percent level. 

Having discussed the system acceptance function, which is 
crucial for all background models, we now return to the princi- 
pal task of a background model: to provide estimates of N « 
and a (Eq. Q & l|2}). Various background estimation tech- 
niques are described below. 

2.2. Ring Background 

A method that is robust in the face of linear gradients in ar- 
bitrary directions is the r/ng-background model. In this model 
a ring around a trial source position (in celestial coordinates) 
is used to provide a background estimate. This is applicable to 
any point in the FoV. The parameter a is approximately the ra- 
tio of the solid angle of the ring (of typical radius 0.5°) to the 
trial source region Q on /O ff, and is typically chosen to be ~l/7. 
However, within the ring the acceptance can not be assumed 
to be constant, since the ring covers areas with different offsets 
from the observation position. Therefore an acceptance correc- 
tion function must be used in the determination of the normal- 



isation a for each position on the ring. The Wng-background 
method is illustrated schematically in Fig.0](left). 

2.3. Reflected- Region Background 

The reflected-region-backgrou nd model was originally devel - 
oped for wobble observations dAharonian et al. H200ll Eo06c). 
but can be applied to any part of the FoV displaced from the 
observation position. For each trial source position a ring of 
"off OFF regions is used (see Fig.|4](right)). Each OFF region 
is the same size and shape as the ON region and has equal off- 
set to the observation position (note that here the ring is cen- 
tred on the observation position, while for the ring background 
technique the ring is centred on the trial source position). The 
method is called reflected-region method because the ON re- 
gion is reflected with respect to the FoV centre to obtain one 
OFF region. In the general case as many reflected OFF re- 
gions as possible are then fit into the ring whilst avoiding the 
area close to the trial position to prevent contamination of the 
background estimate by mis-reconstructed y-rays. Due to the 
equal offset of ON and OFF regions from the pointing direc- 
tion of the system, no radial acceptance correction is required 
with this method and a is just 1 /n Q ff . This is particularly helpful 
for spectral analysis where an energy-dependent radial accep- 
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Fig. 4. Count map of y-ray-like events from 5 hours of H.E.S.S. observations of the active galaxy PKS 2155-304 (Ahar onian et all 
2005d). Note that the data were taken in wobble mode around the target position with alternating offsets of +0.5° in declination. 
The ring- (left) and reflected-region- (right) background models are illustrated schematically. 
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Fig. 5. Illustration of the template-background model. Left: Distribution of the mean reduced scaled width (MRSW) from y-ray 
as compared to proton simulations. This parameter is defined as the deviation in units of standard deviation of the measured image 
width from the simulated expectation value, averaged over all telescopes participating in an eve nt. The s eparatio n potent ial of the 
MRSW is clearly seen, it is frequently used for background suppression in H.E.S.S. analyses dAharonian et alJfeOOSd ). Events 
falling into the Signal region are considered y-ray-like events, events falling into the background region (3.5<x < MRSW < 8cr) 
are considered cosmic -ray-like events and are used for background estimation. Right: System acceptances for the Signal and 
Background regimes as indicated in the legend. The distributions are generated from OFF runs. The Background acceptance is 
normalised to the area of the Signal acceptance in the central 1.5°. 



tance function would otherwise be required. We note that in 
case the y-ray source was observed under a large range of off- 
set angles with respect to the system pointing direction, for ex- 
ample as part of a sky survey, the normalisation a might differ 
substantially from run to run. In this case, a suitable averaging 
procedure has to be applied to both nominator and denominator 
of Eq. 13 the exposure measure is weighted by a factor taking 
account of the offset of the source from the pointing direction 
(this factor might be calculated as the ratio of the y-ray accep- 



tance at the offset of the run to the acceptance at a reference 
offset). 

2.4. Template Background 

The femp/afe-background model was fi rst developed f or the 
HEGRA instrument and is described in Rowell (2003). This 
method uses background events displaced in image-shape pa- 
rameter space rather than in angular space. A subset of events 
failing y-ray selection cuts are taken as indicative of the lo- 
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Fig. 6. Left: Maps of statistical significance for the field around the supernova remnant SN 1006 derived using four different 
background models. 6.3 hours of 4-telescope H.E.S.S. data without significant y-ray signal are used. The black stars in each 
field mark the position of bright (m# = 2.5 and mg = 2.9) stars. The upper left-hand plot is produced applying the field- of -view 
method, for the purpose of comparison, the other plots show the difference between alternative models and the field- of -view 
model. Note that the bins as plotted here are correlated since the significance is calculated integrating events in a circle of 0.1° 
radius of each trial source position. Right: Distributions of significance in the FoV around SN 1006 for the different models. The 
solid black curve illustrates the expected normal Gaussian distribution. As can be seen, deviations from the expected behaviour 
are at the less than 1 % level. 



cal background level. The approach is demonstrated in Fig. [5] 
(left). Events falling into the Background regime are taken as 
OFF counts, y-ray-like events from the Signal regime are ON 
counts. The normalisation a is calculated as the number of 
events in the Signal regime, excluding the source region, di- 
vided by the number of events in the Background regime. A 
correction factor depending on the position in the FoV has to 
be applied to a since the system responds differently to the 
cosmic-ray-like than to the y-ray-like events. Therefore, an ad- 
ditional radial acceptance curve for the Background regime has 
to be determined. This cosmic-ray acceptance curve depends 
on the choice of Background regime. In practice it turns out 
that the system acceptance becomes very different from the y- 
ray acceptance if Signal and Background regime are too far 
apart. This is undesirable because the necessary correction fac- 
tor would vary strongly within a FoV, potentially increasing 
systematic uncertainties. The choice of Background regime is 
thus a compromise between good separation from the Signal 
regime and small a (i.e. reasonable event statistics), and ob- 
taining a background system acceptance function which does 
not differ substantially from the y-ray acceptance. For the par- 
ticular choice of Background regime applied here, the differ- 
ence of the two system acceptance curves can be seen in Fig. [5] 
(right). In the central part, for event directions close to the sys- 
tem pointing direction, the two curves are very similar. For 
offset angles larger than ~ 1.5° pronounced differences occur. 
The ratio of the y-ray and cosmic -ray acceptances, which is 
required to determine a, will not be constant over the FoV. 



The template method has the advantage that the back- 
ground is determined in the same region as the signal and hence 
any localised problem, for example due to a bright star, will af- 
fect both signal and background. Whether the effect is equal 
for both and therefore cancels out depends on the choice of the 
Background regime and has to be checked from case to case. A 
drawback of this method is that exact knowledge, not only of 
the Signal, but also the Background acceptance is required, po- 
tentially increasing systematic uncertainties related to the mod- 
elling of the system acceptance. 

Another method that has been applied to H.E.S.S. data is 
the weighting method which is related to the template back- 
ground. Signal and background are estimated simultaneously 
from the same portion of the sky. Events with directions associ- 
ated with a certain sky bin are assigned two weights, one for the 
assumption that it is a signal event, one that it is a background 
event. Subtraction of the accumulated bin content yields the 
y-ray excess. This approach is not pur sued fur ther, but is d e- 
scribed in detail elsewhere jLemoine-Gourmard et al. 2006). 

2.5. Field-of-View Background 

For the ^zeZc/-o/-v/ew-background model, the entire field (ex- 
cluding regions of known y-ray emission) is used for the nor- 
malisation of an acceptance model to the data and the normal- 
isation a approaches zero. The acceptance model is derived 
from the set of OFF runs mentioned above. Given an observa- 
tion at a certain zenith angle, a model background map is cre- 
ated by rotating the radial acceptance curve (cf. Fig.Q of the 
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corresponding zenith angle band. The advantage of this model 
is that it can be readily applied to extended sources and re- 
sults in the highest possible statistical significance. However, 
the method is sensitive to deviations of the true system accep- 
tance from the model applied. 

2.6. Classical ON/OFF Background 

As mentioned earlier, in traditional ON/OFF mode twice 
the observing time is required for each source, providing a 
strong disincentive for this approach to background modelling. 
Nevertheless, the ON/OFF mode has a powerful advantage in 
that no assumption is made for the system acceptance, except 
that it is the same in both exposures. As the ON and OFF 
runs have identical pointing direction in the horizon system, the 
only assumption made is that the acceptance is not dependent 
on conditions fixed in celestial coordinates, such as stars and 
night-sky background light. This advantage motivates a modi- 
fied form of ON/OFF analysis that has been applied to H.E.S.S. 
data. A sizable fraction of the fields observed with the H.E.S.S. 
instrument contain no significant y-ray signal. These data can 
be used as an archive of OFF data. For a given set of ON runs, 
a set of OFF runs matching in zenith angle is selected from the 
archive. The normalisation, a, between ON and OFF runs is 
deduced from the total event numbers in the two runs, exclud- 
ing the nominal y-ray source region. 

3. Model Comparison 

A satisfactory background model must meet one main cri- 
terion: when applied to many trial positions it should pro- 
duce a normal distribution of significance for an empty field. 
To test this, four of the background models described here 
have been applied to the H.E.S.S. data of SN 1006, which 
(as a lready m entioned) shows no evidence for y-ray emis- 
sion jAharonian et al 12005 eh . Figure [6] (left) compares signifi- 
cance maps of this one field derived using four different back- 
ground models. Hard cuts have been applied and the signifi- 
cance at each trial sky position was calculated integrating ON 
events within a circle of 0.1° radius. A map of absolute signif- 
icance is shown only for the fie Id- of -view model. To ease com- 
parison, for other models the difference in significance to the 
reference field-of-view-model is plotted. The maps show sat- 
isfactory agreement with each other on the 1 cr-level, within 
the expected statistical fluctuations. The difference maps are 
roughly constant throughout the field, with the only excep- 
tion being the regions close to two bright stars in the lower 
right of the FoV. As is described in more detail below, bright 
stars cause a reduction in the local rate of events, producing a 
dip in all significance maps except that derived using the tem- 
plate-background model, in which case the bright stars influ- 
ence also the background estimate. For the particular choice of 
background regime employed here, the template model slightly 
over-corrects for the dip in the y-ray acceptance, producing a 
positive significance at the star positions. 

Figure [6] (right) shows the distributions of significance 
values of each trial source position for the maps shown in 
Fig.|6](left). The regions close to the two bright stars have been 



excluded. The distributions show satisfactory agreement with 
the expected normal Gaussian. We note that deviations on the 
mean significance slightly larger than 1 / VMriais are expected 
due to the correlations between the signal and background esti- 
mates in neighbouring positions. These distributions show that, 
at least for this FoV, the systematic error on the assignment of 
a statistical significance to the signal at a given position is at 
the < 0.1 cr level. Furthermore, the agreement between different 
models suggests that under normal circumstances they all pro- 
vide valid background estimates with fluctuations at or close to 
the expected Poisson behaviour. 

Figure[7]compares a y-ray count map of the field around the 
supernova remnant RX J1713. 7-3946 to three different back- 
ground model map s. The count map is gen erated from H.E.S.S. 
data from 2004 (lAharonian et al.l l2006bl) with four observa- 
tion positions, all offset by 0.7° from the centre of the rem- 
nant. In each observation run, the usable range of the FoV 
is restricted to the central 2° around the observation position 
causing edges in acceptance when overlaying data from differ- 
ent observations. The background models are normalised by 
a and shown on the same scale. It is apparent that the dif- 
ferent models have different levels of statistical fluctuations. 
For the template-background map the statistics are reason- 
ably good, a x 1/14 when choosing a background regime of 
3.5<t < MRSW < 8cr (cf. Fig.|5J and hard cuts, a is practically 
zero in case of the^e/ii-o/-v/evv-background model, statistical 
fluctuations are negligible, consequently the background map 
is very smooth. The OFF-data map (from the ON/OFF anal- 
ysis) has the largest a (» 1), and statistical fluctuations are 
at a considerable level resulting in a comparatively low sta- 
tistical significance of the signal. For illustration contour lines 
are overlaid on all four sky images in the figure. Apart from 
differing event statistics, the three background models are in 
good agreement both in terms of shape and absolute level. They 
clearly provide an appropriate description of the background of 
the y-ray count map shown in Fig.Qa). 

A systematic comparison of the background level estimated 
by different models was performed for the whole H.E.S.S. 2004 
Galactic plane survey. Figure 4 of Aha ronian et alJ J2006a) 
shows the correlation between the background estimate for 
each grid point in the sky derived using ring- and template- 
background models. The correlation is close to linear over a 
large dynamic range. The spread is consistent with statisti- 
cal fluctuations in N g. The slope of the correlation is 1.007 
an d both background estim ates are consistent within 1% - 
see lAharonian et alJ d2006ab for details. 

Another demonstration of the validity of different back- 
ground models is shown in Fig. [8] For two H.E.S.S. data sets 
with significant y-ray signal, slices along Right Ascension and 
encompassing the sources are shown. Overlaid, in both cases, 
are the fie Id-of-view and template background models for these 
data sets. In both cases, at different regions in the sky, for 
an extended and a point-like y-ray source, there is clearly a 
good match between both models and data in regions outside 
the y-ray sources. Moreover, it is evident that the features of 
the gamma-ray morphology of RX J1713. 7-3946 after back- 
ground subtraction are robust and remain unchanged when ap- 
plying different background models. Note that this can also be 
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Fig. 7. Illustration of different background models for the case of the supernova remnant RX J1713. 7-3946. Note that hard cuts 
have been applied here, a) Raw y-ray count map generated from H.E.S.S. 2004 data for the field around the remnant, b) - d) 
Normalised background maps derived using three different approaches (discussed in the main text). Overlaid on all four images 
are white contours for illustration. They are equally spaced at 10, 20, and 30 counts and are deduced from a Gaussian-smoothed 
version of the raw colour images to reduce the impact of statistical fluctuations. 



seen in Fig. 6 of lAharonian et al] J2006bl) which shows a lin- 
ear correlation of gamma-ray excess counts for the sky region 
around RX J 17 13.7-3946 for two different background sub- 
traction methods, the weighting method (mentioned a bove and 
discussed elsewhere JLemoine-Gourmard et al. 12006^ and the 
field-of-view method. 

4. Effect of Stars 

All background estimates presented here rely on the homo- 
geneity of the (y-ray or hadron) acceptance across the FoV. 
While detector acceptance inhomogeneities are typically of the 
order of 3% or less, they may reach higher values in special 
cases such as large zenith-angle observations or in the pres- 
ence of strong sky -brightness variations (most frequently due to 
stars). Figure|5] shows as an example the event rate in arbitrary 
units, as a function of the distance to bright stars in the FoV, for 
different bands of stellar B-band magnitude. The curves were 
derived from the H.E.S.S. Galactic Plane survey dataset, av- 
eraging over all stars in the respective magnitude interval. It 
can be seen that for stars with B-magnitude smaller than 5, 
the event rate at the position of the star decreases noticeably. 



This effect can be explained as a consequence of the automatic 
switching off of pixels when DC illumination reaches poten- 
tially damaging levels. For events in which the shower core is 
located close to one of the telescopes, the event for that tele- 
scope will have a hole in the middle of its Cherenkov image due 
to pixels being switched off by the star. An image with a hole in 
the middle might a) be thrown away by the image cleaning, or 
b) fail the shape cuts. Therefore less events are reconstructed 
in the direction of the star. The effect increases with increas- 
ing brightness of the star. The histograms shown in Fig. [9] are 
derived using std cuts. For hard cuts the effect is less severe. 
This is expected, as the influence of individual pixels that are 
switched off is smaller for the larger images required with the 
increased image amplitude cut of 200 p.e.. Thus, for the hard 
cuts, no significant dip at the position of stars can be seen in 
the magnitude 5 to 6 band. For the magnitude 4 to 5 band, a 
significant dip occurs for both sets of cuts. The magnitude of 
this effect is si milar to that found prev iously for the HEGRA 
telescope array ( Puehlhofer et al. 2003). 
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Fig. 8. Illustration of the agreement between data points and background model (solid and dashed lines). H.E.S.S. data for the 
extended y-ray source RX J 17 1 3.7-3946 (a)) are comp ared to PKS 2155-304, a point-like extra-galactic y-ray source for which 
a high-statistics data set exists ( Aharonian et al. 2005d) (b)). Slices along Right Ascension (RA) through the centre of (and fully 
encompassing) the sources, are shown. Two background models are shown, the fie Id-of-view- (dashed line) and the template- 
background model (solid line). Note that the steps in the distributions are artifacts of the analysis: for each observation run, the 
usable range of the FoV is restricted, and the distributions are produced from different observation positions. 





ring 


reflected- region 


template 


field-of-view 


ON/OFF 


Contemporaneous 


Y 


Y 


Y 


Y 


N 


FoV position 


N 


Y 


Y 


N 


Y 


Sky position 


N 


N 


Y 


N 


N 


Event statistics 


Y 


Y 


Y 


Y 


N 


Event type 


Y 


Y 


N 


Y 


Y 



Table 1. Overview of the properties of the different background models described here. For each feature (described on pagellO> 
we quote if a given model fulfills (Y) or fails (AO this condition. 



5. Choice of Model 

In the following we summarise the advantages and disadvan- 
tages of the available background models and list the analysis 
types they are most suited for. We start by listing properties of 
an ideal background model: 

Contemporaneous: Background events for a given data set 
should stem from a contemporaneous observation period 
to avoid incompatibilities due to ageing effects of the tele- 
scope system. 

FoV position: Background events should be accumulated at 
the same or a similar position in the FoV, meaning that the 
angular distance to the system pointing direction should be 
equal to the signal events (because of the angular depen- 
dence of the system acceptance, cf. Fig.0. 

Sky position: To assure a similar level of night-sky back- 
ground light background events should be recorded from 
the same region of sky. 

Event statistics: To reduce fluctuations, the event statistics of 
the background should be considerably larger than the sig- 
nal one, implying a normalisation factor s« 1. 



Event type: Signal and background events should be of the 
same type, from the same region in image-parameter phase 
space to assure a similar system response and reduce the 
importance of the correct system-acceptance model. 

It is obviously impossible to fulfil all of these requirements 
with a single background model. Any choice can only be a 
compromise. It is therefore important to apply different models 
to a data set thereby cross-checking the results. Table[5]classi- 
fies the models discussed in Section|2]in terms of the properties 
listed above. The advantages and shortcomings of the different 
models are: 

ring background: The model has the advantage of provid- 
ing a conceptually simple prescription for the background 
determination. It is rather insensitive to deviations of the 
actual relative acceptance of the data set from the model 
acceptance function, as it only relies on the relative nor- 
malisation in a limited nearby area around the source bin 
in the FoV. Any linear gradients in the system acceptance 
are averaged out because of the summation on the ring. 
However, when testing larger source extensions, larger ring 
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Fig. 9. Effect of bright stars on the event rate of the system 
for standard and hard cuts. The relative event rate is shown 
as a function of the distance to bright stars for different B- 
magnitude bands. The 200 p.e. size-cut data set has been scaled 
by 0.75 for clarity. 



radii are required, demanding better accuracy for the rel- 
ative acceptance correction across a larger portion of the 
FoV. For the determination of energy spectra, this method 
is disfavoured, since the acceptance curve changes with en- 
ergy (as shown in Fig. [2] (left)). Any attempt to correct for 
this, e.g. by determining acceptance curves in energy bands, 
would introduce another source of systematic uncertainty. 
Additionally, the OFF events have a different distribution 
of offsets from the centre of the FoV than the ON events. 
Since the effective areas used for the spectral analysis de- 
pend on the camera offset, the different offset distribution 
introduces again an additional systematic uncertainty, even 
if one corrects for it. We note furthermore that in the case of 
several y-ray sources or indications for sources in the FoV, 
a case which is becoming more common in Galactic ob- 
servations with the current generation of experiments, the 
ring background suffers from additional systematic uncer- 
tainties, that is, possible y-ray contamination: If the source 
is surrounded by several 3cr spots, some of which might 
be actual so far undetected sources, they would all be in- 
cluded in the background estimate and lead to a systematic 
over-subtraction. 
reflected-region background: It has the advantage that it is 
independent of the exact shape of the acceptance function. 
It simply relies on the assumption of radial symmetry of the 
acceptance. Additionally the distribution of offsets from the 
centre in the ON and OFF events is the same. This makes 
the model especially suited for the background estimation 
for energy spectra. However, this approach relies on a suit- 
able observation strategy, it cannot be applied if the obser- 
vation positions of a data set are within an extended source 
region, or, as mentioned above, in case of too many other 
y-ray sources in the FoV. In this case one either ends up in a 
situation where it is not possible to define a reflected back- 
ground region without overlap with an actual known source 
region, or, in case of close-by sources just below detection 



limit, one might obtain a y-ray contaminated background 
estimate. 

template background: This technique has the advantage that 
it is better suited to largely extended sources (which fill 
a sizable fraction of the experiment's FoV) than the ring- 
background method as long as an acceptance model (say 
from OFF data) is available. It is, however, sensitive to un- 
certainties in the relative acceptance determination between 
the Signal and Background regime across the FoV. Large 
differences of the system-acceptance functions in the two 
regimes potentially increase systematic uncertainties of the 
normalisation factor a. 

For the estimation of energy spectra the background esti- 
mate must consist of events with a similar distribution of es- 
timated energies to the background events of the source re- 
gion. The template model does not meet this criterion since 
the events in the background regime will differ in estimated 
energy from those in the signal regime. Also the problem 
of the acceptance curve changing with energy is present as 
in the case of the ring background but more severely since 
it applies to the signal as well as to the background regime. 

field-of-view background: The model can readily be applied 
to any data set to investigate the source morphology. It is es- 
pecially well suited for very extended sources that cover a 
large fraction of the FoV and yields the maximum possible 
signal significance since the normalisation factor a is prac- 
tically zero. The caveat is that it is sensitive to deviations of 
the model from the true system acceptance. For example, 
pronounced night-sky brightness variations within a single 
field and data sets with unbalanced wobble observation off- 
sets might cause distortions of the true system acceptance. 

ON/OFF background: This classical approach is a robust 
method to perform cross-checks and explore systematic un- 
certainties of spectra of very extended sources. Its advan- 
tage is that it can be applied to any data set, independent of 
the source size and the observation strategy pursued for a 
given source (as opposed to the reflected-region approach, 
which relies on observation positions outside the source re- 
gion). Its caveats are the loss of a factor of two in ON- 
source observation time and possible changes in the night- 
sky background level between ON and OFF data. 

To summarise, the best suited background estimation tech- 
nique for the extraction of various aspects of the y-ray signal 
are: 

□ Source detection: the nng-background model has in general 
fewest systematic biases for this purpose. Only in the case 
of busy sky fields which contain multiple potential gamma- 
ray sources issues due to gamma-ray contamination arise. 

□ Spectral analysis: a reflected-region background is most 
suitable due to the identical offset distribution of ON and 
OFF regions. For very extended sources, however, or for 
extended sources which have not been observed with suf- 
ficiently large wobble offset, or in case there are other y- 
ray sources in the FoV and one cannot define background 
regions, the ON/OFF -background model is the only appro- 
priate one. 
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Fig. 10. The impact of an offset between observation position 
and target source (wobble offset) on the source significance. 
For a reflected-region background an offset is required to pro- 
vide OFF regions at equal offset. The relative significance has 
a broad plateau between 0.5 and 0.7 degrees for point-like 
sources. 

□ Morphology of very extended sources: the field- of -view 
model provides an effective way to investigate the morphol- 
ogy of sources too large to be effectively handled by other 
methods. 

5.1. Optimal Wobble Offset 

Besides choosing a background model best suited for a given 
source and the analysis task, the optimal observation strategy 
needs to be considered. It is governed by the source properties 
and the preferred background model to be applied. To derive a 
background estimate and extract a signal from a single dataset 
it is normally necessary to observe a potential source with an 
offset with respect to the pointing direction of the system. For 
a given background model and for known radial system accep- 
tance, the optimum observation (or wobble) offset opt can be 
calculated. This offset is defined as the one which maximises 
the significance per unit observation time (S / yi) of the source. 
For ring and template backgrounds opt = 0, but this is unde- 
sirable for two reasons. The first is that only if the source is off- 
set from the system pointing direction can the one-dimensional 
radial system acceptance be extracted from the data set under 
study (because of the exclusion of the source region for accep- 
tance generation). The second reason is that a spectral analysis 
becomes more difficult and one introduces systematic difficul- 
ties for extended sources since the reflected-region method can- 
not be applied. In fact, for the reflected-region background, opt 
is a compromise between the number of available OFF regions 
(which increases with increasing offset) and the fall-off of the 
system acceptance for large offsets. In Fig. ^| we explore the 
significance as function of wobble offset. For point sources in 
H.E.S.S., S / yfi exhibits a rather fiat plateau between 0.4° and 
0.7°. For moderately extended sources (<x soul - ce ~ 0.2°) the op- 
timal offset increases to 0.7°- 1.0°. For very large sources the 
best strategy is to observe just outside of the source, with some 
safety margin to account for the finite angular resolution, so 
that one reflected off region is available without y-ray contam- 



ination. Targeted H.E.S.S. observations are typically taken at 
0.5° or 0.7° offsets. 

6. Conclusions 

Several different background models are available for ground 
based Cherenkov astronomy. Different models are appropri- 
ate for different purposes. Searches for weak sources are 
best performed with the robust Wng-background model. For 
spectral analysis the reflected-region background is favoured. 
Extracting the morphology of extended sources is often most 
reliable with the field- oj -view or template models. In general, 
a comparison of several models (with different systematics) is 
necessary to establish the existence of a new source. It is im- 
portant to remember that the estimated statistical significance 
of a source is largely irrelevant if background systematics are 
not under control. 
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